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HIEFACE 



A modern high speed digital computer uses many vacuum tubes in 
bistable circuits. These circuits and hence the computer, are susep- 
tible to the failures of the vacuum tubes. Although considerable 
progress has been made towards improving the reliability of vacuum 
tubes, they do not approach that of the later-developed "solid- 
state" devices. In addition to greater reliability offered by these 
devices, the small size and low power consumption is particularly 
enticing where large numbers of bistable circuits are needed. In 
particular, the development of the ferrites and titanates for high 
permeability and dielectric constants introduced non-linear elements 
in a component size that surpasses the highly lauded transistor. As 
with other new circuit elements much is left to be done, especially 
in improvement of the materials involved. However, the first step 
was taken by C. F. Spitzer fl] in the analysis of the two-branch, 
bistable, series resonant circuit, using a non-linear inductance. 

The next move is the analysis of the non-linear capacitor circuit 
which this paper purports to do. The next obvious step would be the 
Inclusion of both devices in a single circuit, and a short develop- 
ment of this appears in hopes of beginning the solution to this prob- 
lem. 

The analysis is to be done with the thought in mind of obtaining 
a set of design equations for a circuit which might be used as a bi- 
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stable network. Both the two-branch, series resonant circuit, which 
consists of series elements in a parallel network, and the two-mesh, 
parallel resonant circuit, which consists of parallel elements in a 
series network, are analyzed. Consideration of triggering schemes, 
although of great importance, occupies a small portion of this paper 
since it was not the primary purpose to design an actual circuit but 
rather obtain the criteria which might lead to a realizable circuit. 
The production of an actual operating circuit awaits the development 
of materials with lower losses and slightly greater non-linearities. 
The necessary changes and the progress being made on this phase of the 
elements are entirely beyond the scope of this paper. 

The author wishes to express his deep gratitude to Charles A. 
Rosen and Charles F. Spitzer of the General Electric Laboratories for 
their extensive suggestions and guidance in the completion of this 
paper. Without the many "skull” sessions with these gentlemen, this 
paper probably woijld not have been written. 
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CHAPIJR I 



INTRODUCTION 



The circuits to be analyzed in this paper are the two-branch, 
series resonant circuit shown in figure I-laj and the two-mesh, para- 
llel resonant circuit shown in figure I-lb. Both of these circuits 
utilize only one non-linear element. In chapter IV a brief develop- 
ment appears of the two-branch, series resonant circuit using two 
non-linear elements as shown in figure 1-2. 

Inasmuch as these circuits contain non-linear elements, methods 
other than linear network analysis must be employed. There are several 
approaches for analyzing such circuits based upon developments by 
Duffing, Poincari , Kryloff and Bogoliuboff and others 2, 3, h, 
and 87- The choice of one approach over another depends on the accur- 
acy desired, the accuracy obtainable, and the exact problem confront- 
ing the analyst. The method of Dryloff and Bogoliuboff is certainly 
more general than that of Poincarej and both are more accurate than 
that of Duffing. On the other hand, the method of Duffing is far 
easier than the others and yet yields the desired information to a 
sufficient degree of accuracy. All of these methods mentioned are 
analytical and suffer from one inherent disadvantage: they produce 

solutions for only the case under consideration without any thoiaght 
as to the totality of solutions available. A qualitative method, 
known as the topological approach, gives an overall picture of the 
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Figure I-la Toro Branch, Loss-Free Series Resonant Circuit 
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Figure 1-2 A two-branch Loss-Free Series Resonant 

Circuit using Two Non-linearities 
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totality of solutions but yields no specific information on any partic- 
ular case J,, and 8_/ . 

The non-linear capacitors available at this writing had relatively 
large losses, were dependent on temperature, and had small non-lineari- 
ties. Furthermore, there was no available test equipment for measur- 
ing the component parameters at the higher frequency of circuit opera- 
tion. Finally, the determination of the component parameters was based 
on some rather sweeping assumptions and hence accuracy was sacrificed 
ab initio. 

For these reasons, and due to its simplicity, a modified form of 
Buffing's method was chosen. In order to clearly understand the devel- 
opments of chapters II through IV, a brief presentation of this method 
and the reasoning leading to it is given in what follows. 

Consider the motion of a simple pendulum having no losses and not 
under the influence of an external force. The differential equation 
of this system, as shown in figure 1-3 is 




I-l 



The solution is begun by using the substitutions 




Equation I-l then becomes, after separating the variables; 




1-2 



k 



Integrating this expression gives; 

(l^X -h A, 

^ ^ 1-3 

where h is the energy constant of the system. From this equation a 

plot of V versus x is made as shown in figure I-ii. From equation 1-3 

and the plots of figure I-ii, the following observations can be made: 

1) h "^0 or becomes negative, which defines an unreal- 

istic case: 

2) if , the curves are closed, encircle the points 

v=0, x=2nTV ,n being an integer, and hence the pendulum oscillates 
about the angle x r 2ntT j and 

3) if , V is never zero, and thus the pendulum con- 

tinues to rotate about its pivot. 

The heavy plot corresponds to the case K-^ , and it can be 
shown that this case cannot exist. However, this value of h divides 
the totality of solutions: those that are periodic and those that are 
aperiodic . 

This procedure just outlined represents the crux of the topologi- 
cal approach. From it there has been determined, once and for all, 
the totality of solutions in the system, but nothing has been learned 
as to the values associated with any particular case. 

It seems evident that if a driving force were included, then 
whether or not the solutions were periodic or aperiodic would depend 
on the amount of energy this force adds to the original system. 



Obviously, if aperiodic solutions are allowed in the electrical 



Figure 1-3 
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Simple Pendulum 







Figure I-l; 



Velocity versus displacement curves 
for simple pendulum. 
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circuit cases, then, at some time, the currents and/or charges will be 
of sufficient magnitude to destroy the circuit elements. Thus analysis 
can immediately refrain from these solutions and concentrate on those 
in which periodicity is mandatory. 

In order to obtain these peidodic solutions we now consider 
Buffing's method which was first applied to the equation bearing his 
name: 

X -h C ^ '¥CiiX •+ - F Cac.Ujt. 

Note that the function of x appearing in this equation resembles the 
first two terms of the expansion of the sin x function. Hence, it 
might be expected that the solutions are much the same as those of 
equation I-l. Such is found to be the case within the region of 
periodic solutions /~1 and 8_/. 

First, assume no damping, i.e. c = 0, and write the equation in 
the following form: 

X = 1-5 

The iteration now starts with the first approximation (the approxima- 
tion X - 0 being considered trivial), X^=/4^60T, and inserts this 
into the right hand side of equation 1-5 ♦ This obtains: 

C^3U)i, 1-6 

Integrating this expression twice, there results 

the constants of integration having vanished in order to preserve 
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periodicity. 

At this point Duffing made a rather far reaching choice. In lieu 
of the usual procedure of prescribing wand then determining the ampli- 
tude, A]_, of the fundamental, he prescribed z ^ and thus obtained 

LO^- 0-1 + -^^ 3 ^ ~ ‘ 1-8 

It is noted that UJ becomes a single-valued function of A but that 
A could be a triple-valued function of U) . It was for this reason 
that A was prescribed, and hence, U> determined. 

The next step of the iteration procedure would be to substitute 
the solution ^ 

into the right hand side of equation 1-5 and again integrate twice 
letting the constants of integration vanish. This, then, is Buffing's 
method. For rapid convergence a^ and A should be small. 

The modified procedure used in this paper removes the restric- 
tions on a^ and A, but implies that the driving force is small of 
order a^ i.e. we set F - a^FQ. 

This begins by adding to both sides of 1-5 and obtaining; 

X 4W ^ — (U) ^ X -- X Q j ^ Cgc- Uft. t 1-9 

Again, iteration starts ;d.th the approximation and substitutes 

this into the right hand side of equation 1-9. 

This gives: 

- ^<^2 <^3cot. 
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Figure I-5a A, 



Relationship for Doffing ’s 
Equation 




Figinre I-5'o A Response Curves for 

Buffing's Equation 
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Now to perform the next step in the iteration it is noted that the 
coefficient of C^aJt — 0 or else secular terms of the type B"tL cu't' 
arise. 

This procedure gives: 

(u’--a.)A -|aj + = 0 

w'-Q., +|ajA%^ i.ii 

which is precisely the same relation resulting before if is re- 

placed by F. After setting theCcotut term to zero, there results: 

X, + ^ 3 U)t 

which gives as its solution 

X, =• -b (1^3 Cut. 1-12 

As in Buffing’s method the choice is to make - A. The procedure is 
then to substitute 1-12 into the right hand side of equation 1-9 and 
repeat the steps above. The iteration is continued until the desired 
degree of accuracy is obtained. 

It is instructive, and will be useful later, to plot equation 
1-8. Although is the dependent variable in this case, it is plot- 
ted on the abcissa in order to show the resemblance of the response 
curves to the linear case. The plots in figure are for the loss- 
less case, while those in figure 1-6 are for the case where losses are 
considered / 8_7* It might be added parenthetically at this point 
that the plots in these figures bend to the right or left according 
to whether a^ is greater or less than zero respecitvely. If a^ is 
zero, the linear case results and the curves are upright. 
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Figure 

Case 



Figure 1-7 Q-E Relationship for A Non-Linear 

Capacitor 
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If a single series circuit were to be driven with a constant 
frequency source whose amplitude could be made to vary, then operation 
would take place along a vertical line such as that shown in figure 1-5 a. 
It should be noted that as the amplitude of the driving force is slowly 
increased a point is reached where there are two possible values of the 
response, A. Since A represents charge in the electrical case, and 
current is the time rate of change of charge, when this critical 
force is reached, the current is said to jump. If, instead of a 
single series circuit, two such circuits are placed in parallel as in 
figure I-la, and a constant current source used, it would seem plaus- 
ible that when the critical force level is reached, one branch current 
would "jump” up while the other branch cxarrent would "jump" down. Such 
an event does occur. 

The question now arises as to what form f(q), shown in figure I-l, 
shall take. If a charge versus voltage plot is taken on one of these 
non-linear capacitors, a curve such as that shown in figure 1-7 resiilts. 
This particular curve is for one value of driving voltage. In figure 
1-8 a set of such loops is shown for various driving amplitudes. The 
heavy line connects the loop tips. Obviously, no simple relation de- 
fines the Q-E curves of either figure 1-7 or I- 8 • 

Since the real interest is in the value of charge which can cause 
bistability, it would seem evident that the large signal loops are of 
more importance. This is due to the fact that the charges must be 
large enough to cause an excxirsion into the non-linear region at the 
time of "jumping" into bistability. 
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If a single-valued function is desired, it is assumed that the Q-E 
relationship follows the mid-points rather than the upper or lower por- 
tions of the loops. This tends to average the characteristics and has 
the advantage that voltage is zero when the charge is zero. 

The assumptions made this far are the primary limitation on the 
accuracy of any solution. It would seem rather irrational to 
choose an involved functional relationship however accurate it may be. 
Thus the defining relation is chosen as: 

an odd function which is required for the midpoint loci. 

The necessary information has now been obtained in order to 
commence the analysis of the bistability problem. 
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CHAPTER II 



TlVOrBRiUrcH, SERIES RESONAI'IT CIRCUIT 



It is the plan of this chapter to investigate the two-branch, 
series resonant circuit. Investigation first concerns itself with 
the loss-free circuit in order to deterraine the basic relationships 
for bistability. After these conditions have been studied, emphasis 
will shift to the low-loss circuit, in an effort to determine the 
effect of resistance on the bistability relationships. 

SBCnON I: THE LOSS-TREE CIRCUIT 

The equations of the circuit of figure I-la can be written as 
follows : 



Subtracting equation Il-lb from Il-la, there is obtained one 
of the basic equations for the circuit; 




+ e, + 




Il-la 




Il-lb 



L 




4- e,-e, o , 



il-2 
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The other basic equation for the circuit is: 



As in chapter I the defining relations for Cj and 0^^ are 
chosen as 



e, = a, 



Since it was decided that the circuit should have a constant 
current source, the impedance seen by the generator voltage, , 
must be that offered by 0^, otherwise the varying impedance offered by 
the two non-linear branches woiild vitiate the constant current h 3 q>oth- 
esis. Thus is leads Oq and hence, since it was tacitly assumed that 
there were no losses, the first iteration begins vrith: 

^,0 = A ujt II-5a 

sr S j II-5b 

where the second subscript denotes the iteration by wMch the value 
was obtained. 

Equation II-2 is rewritten as follows: 




if- fc- e.) , H-6 

The next step is to substitute equations II-!?a and II-5b into 
the right hand side of equation II-6. This results in (see appendix 
III for expansions of powers of trigonometric f-unctions): 
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II-7 



The conditions for periodicity require that: 

-QAC^B)- 30^ (A^- B^) = <3. 

^ 9-L 

iOisregarding the trivial case of A = B and using equations 
II-3 and II-5 there results two equations relating A and B: 



-<Ps -t A-tB 

= — 5 4 /yr\CA^+AS + 6 ^) ■ 



where 



s= 



/yy\ r 



L > 

- 4S a- 



ii-ii 

11-12 



y CL. ' 



^ is a straight line whereas ^ is an ellipse. Using the rota- 



/ 

tion matrix 



'A^ 




?a 






- 
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( — 



Il-lit 



equation II-IO becomes; 






11-15 



The half axes are 



^<=< " / ) and 11-16 

^/S - y ' 11-17 



For the two-branch series resonant circuit the second iteration 
must be performed in order to obtain bistability conditions . The only 
information obtainable from the above relations, is that in order to 
have real solutions of ^ S > 0 and hence Uj ^ (JJoj where 
means: "must be greater than". This condition can be deduced 

intuitively as follows . For bistability k B and hence one of the 
two identically designed branches must offer less opposition to char*ge 
movement than the other, due to the non-linearity of its branch capaci- 
tance. This requires one branch to be "tuned" nearer to UJ than the 
other. Now for a capacitor defined by an equation such as II -1;, it 
is found that the effective capacitance presented to the fundamental 
voltage decreases as the amplitude of the impressed, fundamental 
voltage increases. Hence as the amplitude of voltage increases^ the 
resonant frequency of the branch must increase. Hence the circuits 
must originally be tuned below the source frequency i.e. CO>UJc- 
Thus, what was intuitively believed has been proven by the inequality, 

s ^ 0. 

Solving equation II-7, after the condition of periodicity has 
been invoked, there is thus obtained 

u>t. ~ II-18 
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The choice as explained in chapter I, is to make A^-A and 
Equating the portions of II-7 identifiable with and q2 respectively, 
there resxilts 

II ^ A >6^ Uj't ~ l^A Soo'C II-19a 

I — II-19b 



where 




a. 



3^Co^L 



/W1 



11-20 



Substituting these relationships into the right hand side of 
equation II-6 there results (see Appendix II for details of the 
mathematical devel opment ) : 






- |03 (^l. 
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Invoking the condition of periodicity and disregarding the trivial 



solution, A = B, there results: 

^ 3 ^ 0 = - S t ^ AB + B^) -H V 

-f J2 jj_22 

Using the rotation matrix of equation Il-lli, equation 11-22 
becomes 

= o = -I, + i (3^ V;«‘) + 

■+^('M'’+35M'l/S''+.3l^y +/S‘). 11-23 



This form is still not too instructive. However, another trans- 
formation reduces the form to where the desired results can be obtained. 
Letting and ^ ^ ^ , 11-23 becomes; 



where 



and 





II- 2 U 




■ -h ^dj 


n-25a 






II-25b 


h(s)^ 


/S' -A (LdC.’Z 6 6 — 





-4- 7 (U<1 

For any value of ^ ^ . has only two finite, real roots (posi- 
tive and negative values of the same magnitude) and hence represents 
a single closed loop.^^ and ^ are shown in figure II -1. The 
developments of these curves and the proof of the single root magni- 
tudes are shown in Appendix II. 

From figure II-l, the mode of operation may be explained. For 
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small values of Eg, and hence , the two series branches behave as 
linear networks and thus A = B. As Es is increased, ^ moves out frora 
the origin parallel to the ^ axis. The operating point remains on the 
intersection of and the line, A - B. In Appendix II it is shown 
that for 0-0 degrees, the slope of the tangent of is parallel 
to the ^ axis. Hence will be tangent to ^ at the operating point 

When ^ becomes tangent to , there are now, two other operating 
points available, P2 and P^. Although no mathematical proof is given 
here, it can be deduced from the work of Stoker ^8*7 that the area 
around P]_ is unstable. Operation, therefore, "jumps’* to either or 
P^. The reasons leading to this deduction of Pq_ being unstable are 
given in chapter III for the first iteration of the two-mesh parallel 
circuit. The reason for postponing this discussion until then is 
apparent from the discussion itself. 

A circuit in which the controlled switching of the operating 
point from P^ to P^,or vica versa^ is obtainable constitutes a trigger 
pair. 

The reader might now infer that the first iteration will not 
predict bistability for the two-branch series resonant circuit. How- 
ever, in chapter III it will be shown that this is an unwarranted con- 
clusionj the same value of Eg which will cause bistability will be 
predicted by two different approaches. 

For operating point P]_ in figure II-l, the following values, 
developed in appendix II, are obtained (O - 0°); 
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5 =^ 




j/ 


11-26 




11-27 


= A+B-:2A 


11-28 


Js> ~ CO (ps 


11-29 



Obviously, if ^ moves out from the origin so far that is fails 
to intersect ^ , bistability fails to exist and both branches revert 
to equal currents. Form figure II-l this appears to occur v^hen is 
tangent to ^ at about 0 - 20°. For this condition, the value of Ig 
becomes 



=r CaJ/0 C>S*' 



11-30 



Greater accuracy can be obtained by finding the value at which 
vanishes, but this hardly seems like a reasonable procedure in 
view of the original assumptions. 

It should be added at this point that after bistability occurs, 
there is no assurance that increasing £3 will increase Ig. This is due 
to the phase relationships possible between A and B. If the two opera- 
ting points are located in the first quadrant of figure II-l, then A 
and B are in phase. If the two operating points are located in the 
second or fourth quadrants, then A and B are I80 degrees out of phase. 
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SECTION 2: THE LOW-LOSS CIRCUIT 



Attention is now given to the case where the two branches contain 
losses as shown in figure II-2. Equation II-6 is rewritten to include 
the effect of resistance as follows: 

-f ( e. -ej . 11-31 

Assuming that R is small compared to Xc^_^ then I3 will still be 
primarily determined by the impedance of C-j_. For the remainder of the 
discussion in this chapter, it will be assumed that A^> B. Hence, I^ 
is very nearly in phase with Ig. Since Ig will lead Eg by 90 degrees, 
the first iteration begins with the approximations: 

II-32a 

~ B (Co"^ II-32b 

Substituting these values into the right hand side of equation 
II-31j there is obtained thereby: 

di- fo + wf 

- ^ P) 

_ -)■ ^ec^Lc^i- i>) 

Lj <Cj 

^ HL 11.33 
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Figure II-2 

Two-Branch Low-Loss Series Resonant Circuit 




Figure II -3 

Phase Relationships for Two-Branch, Low-Loss Series 
Resonant Circuit* 
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In order to presejrve periodicity all of the fimdamental component 



terms of UJ must vanish identically. Hence there results; 

- [jsA4wA^]a^^o2^ -4 4^60 2!:, 

Ehqjanding the terms involving and equating the coefficients 

of sin wt and cos wt respectively, there results; 

ip + 'f - I'SA II-35a 

II-35b 

The phase diagram appears in figure II-3. Since R was presumed 
small, ^ is very nearly zero degrees and hence cos !• Using 

this fact and dividing equation II-35a by I 1-3 5b, there results; 

or 

- s 4r^(A’VAB4B’-)- ” 11-36 

If R is allowed to vanish, the first iteration, loss- free case is 
immediately derivable from the above . Experimentally An is of the 
order 8 x 10”'^ whereas s is unity or greater, so that the assumption^ 
S‘5^W^is not too unreasonable.' If it is allowed that s >^An6^ then 
11-36 reduces to 
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11-37 



Now ^ will not be a real illipse, and hence bistability cannot 



result, unless 



R ^ 



a., 5 

(JO 



11-38 



which places an upper limit on the losses. 
For B ^ 0 11-37 gives 












s/nn 



11-39 



If R is allowed to vanish, the first iteration, loss-free case 
is obtained, as can be seen from equation II -10. 

From figure H-2, it "can be seen that; 

(?o = L +• iP + <2, Ii-lio 

or 

€© = A ii-hi 

The magnitude of Eq can be found by first finding its square as; 

Substituting the value of A from 11-39, and upon simplification, 
there results: 
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Since 




11-hh 






The ratio of A to B (i.e. the ratio of the larger to the smaller 
charge magnitudes) is defined as the switching ratio, ^ , and can be 
obtained from equations 11-39 and 11-kSt 

A — » 

^ - UOR 

Now from equations II-l^ and Il-lib, remembering that B^O, there 
results : 

H-ii7a 

* CL^ B . II-it7b 

Since the usefulness of this circuit is predicated upon a suffic- 
iently large difference between the two branch output voltages, it may 
be desirable to solve for this voltage difference: Hence, by sub- 
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traction, and by using equations 11-39 and there results; 




wMch for the loss-free case gives 



E,--E^= <^i(l # 



II-1;9 



Thus it can be seen that although the switching ratio, , is 
independent of the degree of non-linearity, m, that the voltage diff- 
erence is inversely proportional to its square root. It should be 
noted, however, that ^ is not the ratio of to E^, for the circuit 
under consideration. This value, still independent of m, is given by; 



The smaller the value of m is, the larger the potential difference 



from equation II-li9 and the fact that2l\ is related to the parameter, 
s, by equation II-ij.6. This potential difference, however, is obtain- 
able only by making A larger than B. But since the switching ratio. 



the linear region and A in the non-linear region of the Q-E relation- 
ship. 

However, although ^ is independent of ra, it must be remembered 
that too small a non-linearity results in currents of such magnitude 





11-50 



between the non-linear capacitors is for a given ^ . This follows 



, has been prescribed, both A and B increase, with B remaining in 
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that they may destroy the components before bistability is reached. 

All of the above relations were derived from the first iteration 
of the circuit with small losses. The second iteration is now com- 
menced in order to find what modifications, if any, need to be made 
to the first iteration relations. 

After the periodicity conditions have been invoked, equation 
11-33 becomes 



Again choosing Aq_sA and there results, for the second approxima- 
tions ; 



Substituting these expressions into the right hand side of 11-31, 






IX-5I 



whose solution is 



(9, = (A, -B, )^cot - 





gives: 
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p 




— (A-E) 

h 



i 9 9 4 ^ i * 



> 11-Sh 

where the dots refer to terms involving 3Uftj S’CoZ^'?^'t ^nj^cotr. 
The conditions for periodicity become: 

*" (A”B) •+■ b 3 ^ O 11-55 

H‘Li 






11-56 



From II-55j the phase relationship is determined as 

30 ^ * 

r 3 ^3 

1 

Thus, since 4X4^>^3>^ — 

~ ITSJ CA-E> ' 



11-57 



11-58 



Solving for A-B results in: 







/i-l 


i 


- 


iiS 

CO R 



s ( X 2 8 ^ 

yio /? 

2 % 3 ^ 



I 
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Assuming that equation II-ij.6 does in fact give a reasonable 
approximation forz^ , there results, upon substituting into 11-59 

< — 

* 11-60 



A i 
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Here again since bistability implies a value of ^ greater than 
unity, a small value of m may necessitate a value of B so large that 



the components themselves might be destroyed. 

Substituting equation 11-57 in equating 11-56 and by using the 
first two terms of the expansion offl©a3^, there results, assuming 
^ to be small: 

^ ^ II_6i 

Collecting terms, excluding the trivial solution, A=Bj applying the 
rotating matrix, Il-lit; and using the polar coordinate form, there 

results for 0 ^ 0°: , 

ow r . . / . r.K 

11-62 



^ =. f/.Ae 1 ^— 



Except for large values of R and , l/~ is approximately equal 
toXC as given by equation 11-27. 

Briefly summarizing these results: 

1) in order to have bistability ^ 

where ^ ■¥■ ,8^5 . 

S ■/ / 



2) the switching ratio is given by 

uTk j 



A = 



11-26 

11-27 

1 1 - 1:6 
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3) circuit losses should be such that 



R < 



a , £ 

c-o 



J 11-39 

U) the difference of non-linear capacitor potentials is given by 



E-E = a,(i + s^yX 



•, -a ’ > ' I ^ J 

5) the ratio of the two non-linear capacitor voltages is: 



11-1:9 



|i = a(I+S - Ji); W 



11-50 



6) the non-linearity should be such that 



> >1 






II-60 



does not require B so large that the component ratings are exceeded. 
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CHAPTER III 



TWO-MESH, PARALLEL RESONANT CIRCUIT 

The plan of this chapter is to develop relations for the two- 
mesh, parallel resonant circuit similar to those discussed in chapter 
II for the two-branch series resonant circuit. It will also be shown 
that the requirement of a constant-voltage source in the place of the 
previously used constant- current source gives rise to a different solu- 
tion curve relating, A,B, and Es, the source voltage. 

As in chapter II, the first consideration will be the loss-free 
circuit and then the low-loss circuit. 



SECTION 1: THE LOSS-FREE CIRCDIT 



The equations for the circuit of figure I-lb can be written as 
follows; 



^ ~~ 




<^0 



L 




Ill-la 



Ill-lb 



Subtracting equation Ill-lb from III -la there is thus obtained 
one of the basic equations for the circuit: 
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0 



III-2a 



Since there are no d.c. sources in the circuit, no loss of 
generality occurs if equation III-2a is differentiated with respect 
to time. This differentiation results in 

7 ^ 2^ in-2b 

This equation is the same as was obtained in chapter II for the 
two branch, series resonant circuit. 

The other basic equation of the circuit is: 



O, + 6^ - 6s . III-3 

As with the series circuit, the defining relations for and <5*^ 



are; 






Ill-iia 



(?^ = CL, -h Ill-itb 

In lieu of the constant current source, a constant voltage source 
is desired. Hence C, should be fairly large. Thus, ig is a small 
quantity. For lossless circuits, and both lag ig by 90 degrees 
if ix)'> CUo • The charge will be in phase with the voltage and the 



first approximations can be written as; 

rr A Coo COt^ 






in-5a 



(^20 - B Coo Cut, 111-^b 

Proceeding as in chapter II, the two solution curves relating A 
and B are; 

3h 



= - S -t /m(A^ -hAB -i- s'") , in-6 

%=0 - - + in-7 

It should be noted that the curve is common to both systems, but 
that the straight line, ^ , has now become the cubic ^ . These 
curves are shown in figure III-l, as derived in appendix II. The fact 
that ^ is a cubic as shown, makes it possible to demonstrate bistabil- 
ity from the first iteration. The operating point for low values of 
Eg is on the intersection of ^ and the line, A=B. As is increased 
^ moves out from the^ axis, its asymptotic line remaining parallel 
to the^ axis. When ^becomes tangent to ^ at Pq_, the operating point 
"jumps” from to or Py Mathematical proof that is an unstable 
operating point is beyond the scope of this thesis. As stated in 
chapter II, however, the instability may be inferred from the work of 
Stoker (_ 8_7* The reasons leading to this conclusion of instability 
will now be shown. 

In chapter II it was shoxm that the half axes of the eU-ipse are 






3 J and 






II 1-8 
III-9 






The first of these represents the "distance" fl*am the origin to 
Since ArB at this point: 






/\ - / III-IO 

Substituting this value into equation III-7 results in 
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Es = 4)r15« 



III-ll 



As suggested in chapter II there is another method by which 
equation III-IO and hence equation III-ll may be derived. It is the 
criterion of this method which leads to the conclusion that, is 
unstable. In chapter I it was shown that the frequency for the loss- 
free Duffing equation was: 

(jj^ = <1, -I- * III -12 

Since Buffing's equation is the same as the equation for one mesh 
of the circuit under consideration if F is replaced by the voltage, E, 
across the mesh, the frequency is given by: 

0^^= ^ - -T ' 111-13 

^ ^ L A 

The locus of the points where the curve of Fig. I-5b has vertical 
tangents can be obtained by differentiation of equation III-13 id.th 
respect to A and setting dco/c/A equal to zero. This procedure results 



in; 



(V^ "■% ~ S' ^ ^ ^ 

>-> ^ L 



Ill-lli 

The value of A, for a particular UJ , which lies on the locus 
given by III-ll; is found by solving equations III-13 and III-ll; simul- 
taneously. This results in the value of A, 



A ^ 7^rr) • 



III-IO 

stoker has shown that the region bounded by equation III-ll; is 
unstable / 8_7. It was assumed that both meshes of the two-mesh cir- 
cuit were identically designed, and that the circuit was supplied from 
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a constant voltage source. Thus, when A reaches the value given by 
III-IO, only one of the mesh voltages can "jump” up. The other "jumps" 
down. Hence, is unstable and operation "jumps" to either or P^. 

The value of Eg given by equation III-ll has another interpreta- 
tion which may be obtained by study of figure III-l. If Eg were larger 
than the value given by equation III-ll, then there is the possibility 
of having two sets of bistable operating points, since ^ will inter- 
sect ^ in four points. Thus for no more than two possible stable oper- 
ating points 

Ej ^ M,(i +-§)y^ ' 

This apparent contradiction between equations III-15 and III-ll is 
explained as follows. Bistability results when Es reaches the value 
given by equation III-ll. However, as the amplitude of Is is increased, 
there is no assurance that the value of Es vnll increase due to the 
phase relationships possible between A and B. As stated in chapter 
II, operating points lying in the first quadrant represent the in phase 
conditions, whereas those operating points lying in the second or 
fourth quadrants represent the 180 degree phase conditions. However, if 
Eg can be caused to exceed the value given by equation III -12 after 
bistability has been reached, then four stable points may be possible. 

As Es is further increased, there is a value beyond which bistabil- 
ity will cease and reversion to the A=B condition occurs. This value 
of Es is determined by finding the common point of tangency between ^ 
and ^ besides the one already found at P]_ . For ^and ^ to osculate 
at a given point, the following relations must be satisfied at the point 
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III-16 



~\A ■ ’ 

" ^8 

when Es has been adjusted so that the point of tangency lies on 
then both the necessary and sufficient conditions have been met. 

Applying equations III-I6 and III-I7 to equations III- 6 and 
III-7 there results, respectively: 

i ni-18 

/-h • III -19 

Mvision of these two equations by each other, and simplification 



gives 



(^^AB-j)(A-B) = 6^ 111-20 

The solution A=B is trivial for the present development and hence 



there results; 

^^=0 = A^+^ABi-e"^ - 3 ^ 



III-21 



Application of the rotation matrix 



A" 




'//a 


- 




'<x' 








'/a _ 




jS- 



III-22 



to equation I 11-21 results in 


^ C5"c5<^ ^ 


III-23 


a hyperbola with half-axes. 




III-2U 




III-25 
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The ellipse, ^ and the hyperbola, ^ will intersect if the half 
axes obey the relations 

^ ^ ^ ; or ^ ’ III-26 

From this inequality another criterion for bistability results: 

S j or III-27 

CU® > -^ ai. 111-28 

The two criteria for bistability, i.e. the required values of s 
and Es, thus far derived have been obtained from the first iteration. 

The second iteration should give corrections to these relations. 
Proceeding as in chapter II the two solution curves thus obtained are: 

9‘3S0= 

= - 4rf (AAB) +^b(A^g^) 

-h A/mO^{A^-A8'’ ). III-30 

The curves are shown in figure III-2 Continuing the reasoning of 
the first iteration procedure, bistability occurs when^ becomes tan- 
gent to at P^. Applying the rotation matrix, equation III-22 to 
equation III-29 and equation III-30 gives, respectively: 

9i=0= -i- ± ^ 

III-31 

III-32 



40 



As before, the solution curves are now brought into the polar 
form through the equations ; . This results in: 

^ ^ III-33 






-h'^» 



where, f(0), g(0)^ and h(0) are as in chapter II; and 

^ (d) ^ SCexX($ - 

^ S-tacS - VCa<^^~a 

J^(^) - 



III-35a 
III -35b 
III-35 c 



Figure III-2 shows a plot of these curves. The curve differs 
but little from except for large values ot/») ov ^ . The curv^ 

shown is therefore an approximation. 

The value of Es required for bistability is obtained by substitu- 
ting into equation III-3ii the values: O degrees and 

Where ^ ^ . 

This procedure results in 

Es = 4 ’X 

where ^ is given by equation 111-36, 

There is no easily derived correction for the required value of s. 
Application of equations III-I 6 and III-l? to equations III-29 and 



111-36 

IH-37 

111-38 



h2 



III-30 results in an equation which is very tedious to solve. Con- 
sidering the approximations first made regarding the Q-E relationship 
of the non-linear capacitor, the solution does not seem worthwhile. 
SECTION 2: THE LOW-LOSS CIRCUIT 

As in chapter II, attention is now given to the low-loss circuit 
as shoTm in figure III-3^ with the phase relationships shown in figure 
III-1|. Equation III-2b then is rewritten as: 

Assuming R to be large compared to Xc^^, the first iteration commences 
with the approximation: 

^ ^ Ill-iiOa 

III-i;Ob 

Substituting these into the right hand side of equation III-39, 
and involving the conditions for periodicity there results: 

Using the expansions of cos(ivi-/S) and sin (tot+S) and equating 
coefficients of like functions there results, since cos^f^ !■ 

( a , 6 ^(a-b^) 

Dividing equation III-i[2a by III-i:2b and subsequent simplification 
results in: 



h3 



4 




A 



e 

V 



Figure 111-35 Low-Loss, Two Mesh Parallel Resonant Circuit 



Figure 




Phase Relationship for Low-loss, Two-mesh 
Parallel Resonant Circuit. 









lH-h3 



fltl ^ 1 ^ - 

uf-i:- J ~ s~ ^ (a\ A% - t e'^) 



2.. — 

Since in practice, rmB is of the order of 8 x 10”"^ it can be 
neglected in comparison with unity. Thus equation III-43 reduces to 



- s - 



j (a\- A B- h 3^) 

5- /7T> 



2 2 

Solving this equation for A 4AB-}-B results in: 



a^+ab+b^ - 



S + ^ 



iii-lUt 



111-hB 



Clearly this reduces to the loss-free case as R becomes infinite. 
If the right hand side of equation III~I|.5 becomes negative, the equa- 
tion fails to define a real ellipse. Thus in order to pennit bistabil- 
ity, the following inequality must hold: 



/? > 



CO i 

fT- 



III-h6 



For A is taken from as; 







Solving equation Ill-la for io^ and using the fundamental com- 
ponent only results, after simplification, in the expression; 




III-1+8 






Since B is small, the admittance of the associated mesh is: 



-drLlW 



y- s 






Also, since B is small, there results from equation 

« h 

= A|6 = Vj j or 



B= — 



(SijWfl 






R ( S ) y/ry)(^ - 



III-50a 



I I 1-5 Ob 






Again, defining the switching ratio, , as the ratio of A to B 
i.e. the ratio of the larger to the smaller charge amplitudes, there 
results: 



z\ = 



S-h / 



OuL 



II 1-51 



Now, from equations Ill-lia and 111-lj.b, and setting , there 

results : ^ 

^3 A III-52a 

Ej = III-52b 

Since the usefulness of this circuit is predicated upon a suffici- 
ently large difference between the two mesh voltages, it may be desir- 
able to solve for this voltage difference. Hence, by subtraction and 
use of equations Ill-li? and III-50b, there results; 






h6 



my;-: 






111-514 



FOI* th© O.Q Q— +.VTi c r»fa/^n r*Q o Q J 




# 



Here again, although the switching ratio is independent of the 
degree of non-linearity, />y)^ the voltage difference is inversely pro- 
portional to its square root. 

The ratio of the mesh voltages is given by 



Since the second iteration of the low-loss circuit produces 
algebraic equations which are not readily solvable and the first, 
iteration of the low-loss circuit has become rather involved, it 
hardly seems reasonable to attempt a second iteration in an effort 
to obtain more information. 

In summary, then, the following conclusions apply to the two- 
mesh parallel circuit: 

1) The driving frequency must be such that for bistability to be 
possible, 

tu’ > OJ,'-. 111-28 

2) The minimum value of voltage required for bistability is, 
from the first iteration. 





3) The losses must be such that 



R > — ; 




III-I 46 



hi 



h) 



The switching ratio is 





S4- 1 



1 




5) The voltage difference is given by 











6) The value of Eg required for bistability is, from 



iteration 






where 




III-51 
I II -53 

the second 

111-38 

111-36 



and 



III-37 










i 



I 

i 

I 
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CHAPTER IV 



USE OF TWO NON-LINEAR ELMENTS 
FIRST ITERATION 



In chapters II and III, the two basic trigger-pair circuits 
utilizing a non-linear capacitor have been briefly studied. The two- 
mesh series resonant circuit, using a non-linear inductance, has been 
analyzed by Spitzer /~ 7 _ 7 "* Intuitively it is felt that the inclusion 
of both types of non-linearities would reduce the currents required 
for the same switching ratio. This chapter will investigate this 
possibility by analyzing the loss-free, two-branch series resonant 
circuit using two non-linearities as shown in figure 1-2. Due to the 
complexity of obtaining a solution and the inherent inaccuracies pro- 
duced by the assumed current and charge functions, the analysis will 
be confined to the first iteration. 

For reasons given in chapter I, the current- flux and charge- 
voltage relationships, are chosen as; 




IV-1 



Furthermore, since 




IV-3 



equation IV-2 can be written as 



h9 



IV-ii 



iv-5 



IV-6 



One equation for the system is given by 

/\! ~^(e, ^ oj 

whereas the other equation is: 

, f 

'^/ -^2 ~ 

The branch currents and voltages are from equations IV-1 and 
IV-2 respectively:* 

-< - ^ -^3 5 ^ 

'S’, = 

^ = a,^ 7^ a^y 

Again, the capacitor, Cij is to be the primary impedance offered 
to the driving voltage. Hence Is leads Es by 90 degrees. For the loss- 
free circuit, the branch currents will be in phase with the source 
current. Is* Thus the branch fluxes will also be in phase with Is^ 

and the first approximation is taken as: 

CoZl IV- 9a 

<y>xo ~ ^ IV- 9b 

Equation IV-5 can be differentiated with respect to time without 
affecting the solutions since there are no d.c. sources. After per- 
forming this differentiation, the equation is solved for — <5^ )• 



IV-7a 

IV-7b 

IV-8a 

IV-8b 
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In the resulting equation, is added to both sides, and there 

is thus obtained; 

~ y^(e, -ej). iv-io 

The iteration begins by substituting the approximations, <?^oand 9 ^ 
into the right hand side of equation IV-10 and then invoking the condi- 
tions for periodicity. 

The most expedient way to do this is to write <2. as 



e -- a, ^ ^ • 



IV-11 



The value of ^ is then obtained from equation IV-3 and the result 
cubed. The value of is then differentiated directly. This proce- 
dure eliminates the tedious progress of computing and then m\ilti- 
plying by ^ . 

From equation IV-3, using equations IV-9a and IV-9b there results; 

^ xj-i— ^•'2' 



Letting 



and 






IV-13a 



IV-13b 



there results upon cubing IV-12; 
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-4^ 









9^7!-, 



IV-lli 



The condition for periodicity hence becomes 






IV-15 

The terms F and G are replaced by equations IV-13a and IV-13b 



respecitvely and upon simplification there results: 



-|S4,^Y/)l6V 



gl 



where 







wf =■ 



. I -w I 



IV-16 

IV-17a 

IV-17b 

IV-17C 

IV-17d 

IV-17e 

IV-17f 

IV-17g 



i-.C // 

The result of substituting equations IV- 9a and IV- 9b into IV-6 is: 



IV-18 



^2 







K 









This is the equation for ^ found in chapter III. Omitting the 
trivial solution, A ^ B, equation IV-16 defines a solution curve, ^ . 
Bistability results when ^ becomes tangent to ^ . This occurs at a 
point on the A = B line. To find the value of Ig required for bi- 
stability, equation IV-16 is first transformed into the polar coordin- 
ate form by the same transformations used in chapters II and III. The 
resulting equation is solved for^^, , and this value is used to compute 
A. Equation IV-18 then produces the required Is. 

The transformations indicated above reduce equation IV-16 to 

where j are defined as before, and 

^ -iVCficfe y-/, IV-20 

o 

If the term,^, is neglected, equation IV-19 is of the same form 
as the equation defining ^ . Since the functions involving the co- 
ordinate $ , are always positive, the effect of the second non-linear- 
ity is to reduce the required current for bistability. The amount of 
reduction depends on the value of l,m, and n. 

There is no ejqjlicit statement derived, for the value which s must 

take. However, from previous developments of the single non-linearity 

M 

circuit, it can be assumed thatS>0, at least. 

Substituting the value 0=0 degrees into equation IV-19 results 
in: 
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X 




V- - 6, IV-20 

This equation is a quartic in ^ ^ . The general form of the 
quartic is written as 

X^-1-a.x^ - i-JtcSi- C. 1C 7^J= a. iv.^ 

where a,bjC and d are real^ constant coefficients. 

To find the roots of equation IV-21, it is necessary to find any 
real root of the cubic: 

c/(a = (7. IV-22 

The roots of equation IV-21 are then taken from the two quad- 
ratics: 

Ef + {y^)~c iv-23a 

xv Lf -4^yc iv-23b 

where ^ is the root found in equation IV-22. 

Obviously this is a tedious process. Fortunately, it needs to 
be done for only one value of 0. However, to arrive at a general 
solution in terms of the parameters l,m,n, and s is not within the 
realm of this thesis. For any given components and source frequency, 
the value of ^ obtained through the quadric solution is substituted 
into the polar form of equation IV-17: 






17-2h 



This procedure gives the required value of Ig for bistability. 

In summary, the results of analysis of the doubly non-linear 
series-resonant trigger pair are as follows; 

1) The value of Is should be less for bistability in the double 
non-linearity circuit than for the single non-linearity circuit; 

2) There is no explicit development for the value which s should 
take, but it is implied by previous study that S 0. 

3) The required value of Is is obtained from auxiliary quartic, 
cubic and quadratic equations, as given by equations IV-20, IV-21 and 
IV-22 . 
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CHAPTER V 



EXPERIMENTAL RESULTS 

The experimental results obtained with the circuits under consid- 
eration are given in this chapter. They are not offered especially to 
confirm theoretical developments, but to demonstrate that the circuits 
are operable. If the experiments agree with the theoretical develop- 
ments, so much the better. 

The charge -voltage relationships were determined through the use 
of the test equipment arrangement shown in figvire V-1 for 60 cycle 
supplyj and by the arrangement shown in figure V-2 for the $ to 20 
kilocycle range. Sample hys teres igrams for frequencies of 60, 5000, 
10,000 and 20,000 cycles are shown in figure V-3. In figure V-li is 
the hysteresigram taken at 60 cycles, but with the vertical gain and 
intensity set so that the loop characteristics are distinguishable. 

The locus of the tips of a nest of 20,000 cycle major hysteresis 
loops is shown in figure V-5. 

For all experimental work the MUCCN type VSE non-linear capaci- 
tors were used. The values of the coefficients for these type capaci- 
tors were obtained from figure V-5 by approximating the average initial 
and final slopes. The Q-E relationship resulting from this procedure 
is shown as the dashed curve in figure V-5. 

The coefficient values thus obtained are 

a , - 9x^0 ^ volts/ c oul V-la 
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Figure V-1 Test Equipment for obtaining hysteresigram at 
60 cycles. The VTVM should be peak reading. 




Figure V-2 Test equipnent for obtaining hysteresigrams from 
5 kcs to 20 kcs. The VTVM should be peak reading. 
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f: 20 kcs 

Emax: 350 volts 

Qmax: 0.27 uco\il 



f: 10 kcs 

Emax: 350 volts 

Qmax: 0.3 ucoul 




f: 5 kcs 

Emax ; 300 volts 

Qmax: 0.299 ucoul. 



f : 60 cycles 

Emax; 375 volts 
Qmax: 0.315 ucoul. 

(Q axis reversed) 



Figure V— 3. Hysteresipr^ms of MITCON type VSE non-linear capacitors 




Figure V-4. Hysteresigram of MU CON type VSE non-line>^r capacitor 
showing open loop for a source frequency of 60 cycles. 
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V-lb 



/W - 7.63^ /O'^ 

V-lc 

The circuit of figure I-la was first tested using the following 
component values: 



raST NO. 


1 


2 


3 


L 


110 mh 


1.62 mh 


0.1(6 mh 


F 


20 kcs 


233 kcs 


233 ks 



For these components, the computed and experimental results are 
tab\ilated: 



rEST 


1 


2 


3 


Cot 


8.l8riO^ 


5.55x10^^ 


1.96x10^2 


S 


0.9h 


2.86 


HI 

On 

O 

• 


Ig (comp) 


73 itia 


1.59a 


I(81( ma 


Is (exp) 


70 ma 


.51^4a 


327 ma 



Since the parameters, a| and a^, were not separately determined 
for the higher frequency, it is not suprising that the experimental 
and theoretical resiilts disagree. However, from tests two and three, 
it was noted that the higher value of s gave a quicker '*o\omp" but that 
the temperature of the crystal oven in which the capacitors were 
located rose 10 degrees centigrade suggesting a radical change of the 
capacitor parameters as a function of temperature. 
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Figure V-6. Ej and E 2 for the two-branch series resonant circuit. 
The source voltage is larger for the waveforms on the right. 





Figure V*7. £1 and E 2 for the two-mesh parallel resonant circuit. 

The source current, I,, is larger for the waveforms on the right. 




Figure V-8. Relative switching ratios end switching rates of the 
single non-linear resonant circuit and the double non-linear cir- 
cuit. Figure e shows the switch for the two-branch series resonant 
circuit using a non-linear inductance! figure b shows the switch for 
the two-branch serlea resonant circuit using both non-linear induct- 
ance and non-linear oapaoitance. 
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Figure V-9. Output voltage for various trigger rates of one branch 
of the two-branch aeries resonant circuit utilizing two non-linear 
elements. The voltage waveform is taken across one of the non-linear 
capacitors. The trigger pulse was 6.1 microseconds in duration and 
65 volts in amplitude. 
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Switching was not obtainable by use of a pulse generator, but 
only by the use of a long-time (five seconds or more) short circuit 
across the high current capacitor. This indicates that the capacitors 
were heating excessively and that the subsequent change in parameter 
values was detrimental to switching. Methods of maintaining both 
capacitors at the same constant temperature were not feasible because 
of limitation of time available. Figure V-6 shows the output voltages 
and E2 at 233 kcs after the bistability condition has been obtained. 
The phase relation conditions should be noted. The in-phase condition 
was first obtained, and by increasing Es, the out-of -phase condition 
resulted. 

The parallel circuit was next investigated at the frequency of 
20 kilocycles. The data on this test are given below. 



TEST NO. 


k 


L 


110 mh 




9 . 18 x 10 ^ 


s 


0 . 9 ii 


Eg (comp) 


ii 82 volts 


Eg (exp) 


578 volts 



Figure V -7 shows the output voltages, E, and E2 after the bi- 
stability condition has been obtained- Note that the two voltages 
begin out-of -phase and approach the in-phase condition as Is is 
increased. 

IkOsses in both the series resonant circuit and parallel resonant 
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circuit cause the large deviation from the 180 degree relationship. 

The fifth test concerned itself with ».he double non-linear, two- 
branch series resonant circuit versus a single non-linear two-branch 
series resonant circuit. The nortrlinear inductors were those used by 
Spitzer 7 • The parameters of these were: 



^ ^ 2 X 10^ araps/oe/) 


V-2a 


^=r 2 X 10^^ amps/ocf>^ 


V-2b 


/T7 = 7.5 X 10^^ 


V-2c 


N = 200 


V-2d 



The capacitor replacing the non-linear capacitance for the 
single non-linearity circuit was chosen to give approximately the 
same switching ratio. 



The results are given below; 



Test 5 


Nonlinear Capacitor 
Linear Inductor 


Nonlinear Capacitor 
Nonlinear Inductor 


F 


233 kcs 


233 kcs 





h92 raa 


li?5 ma 



Figures V-8a and V-8b show the relative switching ratios and 



switching rates of the two circuits. The double non-linear circuit is 
shown to switch in about two cycles of the supply voltage whereas the 
circuit using only the non-linear inductance switches in about two and 
one-half cycles. 

The sixth test consisted of operation of the double non-linear 
two-branch series resonant circuit at various switching rates. The 
results of this are shown in figure V-9« A Hewlitt-Packard Model 212A 
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Pulse generator was used as the trigger source. 

Only during a brief interval of time as switching obtained on 
the two-branch series circuit using a non-linear capacitor. It was 
apparent from the experiments that the temperature change of the capac- 
itors is the primary deterrent to switching. 

It is felt that the data presented confirms the theory of the 
first four chapters to a good degree. The design equations seem to be 
valid for operation near the frequency region in which the capacitor 
parameters were determined. At high frequencies, the losses are so 
great in comparison to the degree of non-linearity, that the circuit 
does not operate within the theory propounded by the analysis. 
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APPENDIX I 



TRIGCaaiING SCHEMES 



It was stated in chapters II and III that a circuit, in which 
controlled switching is obtainable between the two bistable points, 
would constitute a trigger pair. Obviously the circuit must also 
include means for obtaining this triggering. The purpose of this 
appendix is to present a few ideas as to how this may be obtained. 

The discussion will be limited to the suggested methods without too 
much thought as to their output requirements. One important question 
is left unanswered viz. by what method is memory produced. This is 
an important question, but is considered beyond the scope of this paper. 

When the pendulum was studied in chapter I, it was found that the 
division between periodic and aperiodic solutions was dependent on the 
energy of the system. The circuits studied in chapters I, II, and III 
were shown to pass from the unistable condition to the bistable condi- 
tion when a particular value of Es or Is was obtained. Although the 
implication was not forthright the approximation of the Q-E functions 
of these circuits to the sine function used in the pendulum case leads 
to the conclusion that the two states-unistable and bistable -are con- 
cerned with the amount of energy possessed by the individual branches 
or meshes. 

As an approximation, only the energy contributed by the fxinda- 
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mental ccsmponent will be considered. 

Consider first, the two-branch series resonant circuit- There 
are two possible methods of trigger application: a pulse applied to 
the capacitor; or a pulse induced in the coil. In figure AI-1, both 
methods are shown. 

From equation 11-37, the value of A for which bistability occurs 
is given by 



A -L 

/i= 7^ f /rn aI^~S/rr} Ai-l 

and from equation 11-39 the value of A for the maximum A in the bi- 
stable condition is given by: 




(Z,^SiW 



AI-2 



The current in one branch is<J<^^and for each of the two states- 
unistable and bistable - it is given by: 

J, = «-3a 

Since the energy is associated with the magnetic field for part 
of the cycle of the driving frequency, the energy contained by the 
branch can be written: 

for the unistable state, or 
for the bistable state. 



AI-U 

Al-i:b 
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Figure AI-1 Two Trigger Schemes For Two Branch 
Series Resonant Circuit. 
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AI-5a 



The energy difference is thus given by: 



L|yv\<L^'“ U “ 



3^ a,^s J 



by using equations AI-1 and AI-2. This is the amount of energy that 
this branch must lose before it is possible to switch it into the lower 
charge state. 

Since it was assumed that B '^6 , the amount of energy this cir- 
cuit must gain by virtue of a pulse is given from equation Al-l^a as 




When branch "1" has lost the amount of energy specified by equa- 
tion Al-^a and branch ”2" has simultaneously gained the amount speci- 
fied by equation AI-6a, switching can occur. In order to assure 
switching, the equality should be replaced with the imequality sign 
as follows 



LW-U > ‘^ (S - a?i‘) 

D > (s - 

Since the trigger pulse is applied to both circuits, it must 
cause this exchange of energy to be as specified. At this point the . 
investigation should inquire about memory. Although the subject 
memory is beyond the scope of this paper, a short attempt will be made 
to explain how the energy exchange may be caused. 

The trigger pxilse adds energy to both branches, *'l” and “2”. The 
addition of energy to branch ”1” does not change the effective capaci- 



70 



tance of the non-linear element by much since it was already operating 
in the saturated region. However, the energy addition to branch "2" 
does cause a large change in the effecting capacitance of the non-lin- 
ear element. Hence the opposition to the fundamental component de- 
creases to the point of non-linear resonance in branch '*2” but does 
not decrease appreciably in branch "I”. Since the source is provid- 
ing a constant current, the current now begins to divide between 
branches ”1” and ''2”, adding to branch ”2” by taking from branch "1”. 
Branch ”2” successively offers less opposition to the current whereas 
branch ”1" successively increases its opposition. The process cascades 
until switching has occurred. VJhen the energy in branch "1" has 
reached the value specified by equation Al-l^a, the trigger pulse should 
cease. Otherwise it would be possible to have branch ”1" return to its 
original state of heavy conduction. The reason for this is that the 
pulse is supplying a certain fixed amount of energy which must divide 
between the two branches. Since the effective capacity of the unsat- 
urated branch is greater than that of the saturated branch, the unsat- 
urated branch is more receptive to the energy. However, when the 
branches reach the energy level specified 'ty equating AI-iia, the recep- 
tion is equally well for both branches. Only the effect of cascading 
can then cause the switch. 

The least amount of energy that the pulse must give is that speci- 
fied by equation AI-6a. Since the energy will divide, andat some time, 
equally, perhaps the pulse energy required should be: 
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AI-6c 



Up= > 



i^L 

3/>v\ 







)■ 



The required energy expression for the two-mesh parallel reso- 



nant circuit can be similarly obtained. 
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APPENDIX II 



DERIVATION OF THE SECOND ITERATION SOLUTION CURVES 



The purpose of this appendix is to show hovr the second iteration 
solution curves were obtained. The second iteration begins with the 
approximations ; 

3urt All-la 

^ 3Ufi . All-lb 

These are substituted into the right hand side of equation II-6, which 
after omission of the trivial solution, A=B, results in: 



First, the linear transformation; 



'A' 




"'/E - ‘/e 








r 

1 





is applied. This results in 



AII-3 



A ^( 7 ‘>< ‘°-tSs^‘je'+AU)3 A 

This equation is next transformed into polar coordinate form 
by the relations and ^ • This results in 

(see appendix III for expansion formulas); 
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AII-5 



where 



Si-h Cd<i^6 

^(^) = 3.Si- a (t^^6 -0.6‘(Uc^^ 

/gcacm^ 



AII-6a 

AII-6b 

AII-6c 



Equation AII-5 is now reduced to a cubic through the relation 
psfl> This results in: 

0= AII-7 

From equations AII-6a, AII-6b, and AII-6c, it is noted that, 
for all 6, > 0^^(a)>0 and > o . Thus the terras of 

equation AII-7 varying in sign only as does fl» 

The coefficient signs are 

Mr^y. — — d- — > 

Thus there is one or no possible positive roots j and there ai‘e 
two or no possible negative roots. Since there are three roots, it 
can be stated that there are two negative roots and one positive root. 
The interest lies entirely with the positive root since the negative 
roots produce complex values of ^ which have no meaning as far as the 
circuit is concerned. The "value" /)=»<b does not satisfy the equation. 
Therefore ^ is a single closed curve. 

It now becomes necessary to find the real root of 



1 



I 



AII-0 



The general form of a cubic equation is 

a, x'-+ sc X a. 

The trigonometric method of solving cubics is used and the 
real root is given by 

where 

[-i 

if O and ^ 

if 2^0 and 

Ml 

and ^ — (ZC — ^ ^ 

^ f30.C^ -a^c/)-y^^ 

The upper sign is chosen in bhe^, evaluation ifA>C and the 
lower sign if h^O- It turns out that for all ^,f>0 and thus 
equation All- 9a is the vehicle by which the roots of ^ are obtained. 
The equations of are tabulated below for values of 0 . 



0 








A'3 


(O 




0° 


-s/m 


■| 1*5 r 


41.26'd 




+1.75 


I)2r3 


10° 


-s/m 


■f l.ii7r 


+1.0 D 


r2 


+1.86 


2 3 
D r^ 


20° 


-s/m 


■fl.385r 


+1.23D 


r2 


+2.065D^r^ 


O 

o 


-s/m 


+1.25 r 


+1.13 ] 


3r2 


+2.22 


D2r3 



AII-9a 

AII-9b 
A II- 9c 
AII-9d 

AII-9e 

AII-9f 
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4^ 




Q 



% (r ) 



HP 


-s/m 


+1.09 r 


+ 


1.07 D 




2.15 dV 


Hf 


-s/iti 


f ^ 


4 


1.00 D 




2 r3 


50° 


-s/m 


+0.913r 


+ 


0.905' D 




2.07 D^3 


60° 


-s/m 


40.75 r 


+ 


0.687 D 




1.28dV 


O 

o 


-s/m 


+0.6l7r 4 


Q.H-70 D 




0.753D^r^ 


o 

o 

OD 


-s/m 


+0.53 r 


4 


0.3ID 


■ 4 


0.37 


90° 


-s/m 


+0.5 r 


+ 


0.25'd tP- 


■ 4 


0.25 



From these equations, the value of^ is obtained for each angle 
specified 





yO ) 




UJ 








,9/^s -h.Sos- 
S +! 






J83Z+.S83 

S-^z 






1* 


. 9jTS +S25 
S-P/ 



± ^ j/JS' .alyji M.-.197] 



^2 . S’?/ s -^.7S. 

5tV 



76 



e 
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•> ^ 




These are the only values needed to form a plot, since ^ is 
symmetrical about both the exi axis axis as can be seen from 

equation AII-U* 

The plots in figures II -1 and III-2 were obtained from the 
above tabulations for an s value of 0.68. Although not labeled, 
the coordinates are in units of m 

The development of ^ and ^ is done in a similar fashion. 
Tabulated equations of , for values of 0 are given below. Since ^ 
differs from ^ but little, except for large values of ^ or m, it 
has not been tabulated. The only real interest lies in the fact 
that both ^ and ^ have only one real positive root. For the value 
0-90 degrees, all coefficients of ^ vanish. Hence, /p andjj^ are 
open curves. Since they were derived from monotonically decreasing 
functions, ^i(^) tiy a linear substitution, they will 

also be monotonically decreasing functions about the axis of symmetry. 
Hence the shape is as shown in figures III-l and III-2. 
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Using the method previously described for finding roots of the 



cubic, there arise solutions such as 




AII-10 
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I 



"fi 




where b and c are constants determined for each value of 0. For 

0 A — B. 

values of ^ ^ ^ , the approximation 

ySvT " ) ^ ^ ah-ii 

can be used. Since the approximation specifies | , there is the 
series equality*. 









A 11-12 



Using only the first term of AII-12, equation AII-10 becomes 

X, - 



An-13 



Again since , the hyperbolic sine function can 

be replaced by the first term of its power series expansion. Thus 



^ sa. 



’I 



All-lii 



For the values of m and a^ used in experiments, Es should be 
less than 100 volts for all the approximations to be valid. It was 
on this basis that ^ and ^ were plotted in figures III-l and 
m-2 respectively. It is not anticipated that a marked difference 
would result if the value of Es were such that equation AII-10 must 
be used. 

The slope of the ^ curve is obtained by application of the 
formula 



- 






A 11-15 
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where T is the angle the tangent makes with the horizontal axis. 
Since the horizontal axis is the odaxis, for the tangent to be 



parallel to the ^ axis, the denominator of AII-l^ must vanish, or 
the numerator must become infinite. 

Differentiating AII-5 with respect to 0, there results; 

- [ ^ Aii-16 

From equations AII-6a, AII-6b and AII-6c, it can be seen that 

g'(0) and h'(0) each have the factor sin 0, sin 20, or sin 

/ 

UO. Thus,^^ vanishes for 0-0 degrees and therefore so does the 
denominator of AII-15* It has then been proven that the tangent of 
is parallel to the axis. 
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APPENDIX III 



MATHEMATICAL FORMULAS USED IN ANALYSIS 



This appendix lists, for the readers' reference, mathematical 
formulas used in the analytical portion of this paper. 

The trigonometric expansions are obtained by applying the bi- 
nominal expansion to the exponential form of the sine and cosine func- 
tions. Thus there results 




Alll-la 
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AIII-U 



In determining that ^was an ellipse and investigating if it 
became an imaginary ellipse due to the presence of losses, the follow- 
ing properties of conic sections were used. 

The general form of a conic is written as 

AIII-5 

The relations, D, d and c are defined as: 



D = 



a h g 
h b f 
g f c 



AIII-6a 



d = 



a h 
h b 



AIII-6b 
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a 4 b 



AIII-6c 



The conic is: 

1. a parabola if d ^ 0, O3 

2. parallel lines (may be coincident or imaginary) if 

d=0, D=03 

3. an ellipse if d>0, eD<0j 

i;. no locus (imaginary ellipse) if d > 0, eD>0; 

5. a point ellipse if d>0, D=r0j 

6. an hyperbola if d<" 0, 0 

7. two intersecting lines if d^O, D=' 0. 

In order to remove the xy term of the generalized conic equation, 
the rotation matrix 




AIII-7 



is applied, where 0 is given by 
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